NASA Technical Memorandum 106789 




Modeling of Stress/Strain Behavior of Fiber- 
ReinforcedCeramic Matrix Composites 
Including Stress Redistribution 



The University of Toledo 
Toledo, Ohio 



Pappu L.N. Murthy and Christos C. Chamis 
Lewis Research Center 
Cleveland, Ohio 


December 1994 



National Aeronautics and 
Space Administration 


(NAS A- TM- 106789) MODELING OF N95-17659 

STRESS/STRAIN BEHAVIOR OF 
FIBER-REINFORCED CERAMIC MATRIX 

COMPOSITES INCLUOING STRESS Unclas 

REDISTRIBUTION (NASA. Lewis 
Research Center) 11 p 

G3/2A 0035005 




MODELING OF STRESS/STRAIN BEHAVIOR OF FIBER-REINFORCED 
CERAMIC MATRIX COMPOSITES INCLUDING STRESS REDISTRIBUTION 


Subodh K. Mital 
The University of Toledo 
Toledo, Ohio 43606 

Pappu L.N. Murthy and Christos C. Chaims 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

A computational simulation procedure is presented for nonlinear analyses which incorporates microstress redis- 
tribution due to progressive fracture in ceramic matrix composites. This procedure facilitates an accurate simulation of the 
stress-strain behavior of ceramic matrix composites up to failure. The nonlinearity in the material behavior is accounted for 
at the constituent (fiber/matrix/interphase) level. This computational procedure is a part of recent upgrades to CEMCAN 
(CEramic Matrix Composite ANalyzer) computer code. The fiber substructuring technique in CEMCAN is used to monitor 
the damage initiation and progression as the load increases. The room-temperature tensile stress-strain curves for SiC fiber 
reinforced reaction-bonded silicon nitride (RBSN) matrix unidirectional and angle-ply laminates are simulated and 
compared with experimentally observed stress-strain behavior. Comparison between the predicted stress/strain behavior and 
experimental stress/strain curves is good. Collectively the results demonstrate that CEMCAN computer code provides the 
user with an effective computational tool to simulate the behavior of ceramic matrix composites. 


INTRODUCTION 

Ceramic matrix composites (CMC’s) are currently a subject of a great deal of research interest due to the 
requirement for new high temperature materials for use in high speed engine structural components and other applications. 
Continuous-fiber reinforced ceramic matrix composites offer promise in this area and hence considerable effort is being 
devoted for their development They offer several advantages such as high specific stiffness and strength, higher toughness 
and nonbrittle failure as compared to monolithic ceramics, environmental stability and wear resistance for both room and 
elevated temperature applications. 

There are several differences between ceramic matrix and, say, polymer matrix composites. In polymer matrix 
composites, the fibers provide stiffness, while the polymer matrix serves primarily as a binder and generally has a very low 
modulus. The ratio of Young’s modulus of the fibers and matrix is typically of the order of 0(10 ) in these composites. 

In addition, to achieve high stiffness it is usually desirable to have a strong interfacial bond between the fiber and matrix 
materials. On the other hand, a major role of fibers in ceramic matrix composites in addition to providing higher stiffness, 
is to provide mechanisms to enhance toughness as the matrix material is quite brittle and fails at relatively low strain levels. 
The ratio of the Young’s modulus of the fibers and matrix materials in ceramic matrix composites is typically of the order 
of 0(10°). It is desirable to have a weak interfacial bonding between the fibers and matrix in CMC’s to enable such 
toughening mechanisms as fiber debonding, fiber bridging, fiber pullout, crack deflection etc. in order to eliminate cata- 
strophic failures and to provide some ductility. Therefore, in spite of some similarities with polymer matrix, the analysis of 
ceramic matrix composites requires an entirely different approach. CEr amic Matrix Composite ANalyzer (CEMCAN) has 
been under continued development under the sponsorship of Hl l bMP (Hig h Temperature Engine Materials Program) at 
NASA Lewis Research Center in Cleveland, Ohio to achieve this goal. CEMCAN is based upon composite micro-, meso- 
and macromechanics and has been developed specifically to simulate aspects unique to CMC’s. It incorporates a novel and 
unique fiber substructuring technique that has many advantages including a more accurate micromechanical representation 
of interfacial conditions, both around the fiber and through-the-thickness of the laminate, and provides a greater detail in 
stress distribution within a ply at a very high computational efficiency as compared to numerical analysis techniques used 
for micromechanical analyses. Recently, a nonlinear analysis capability which includes the local stress redistribution due to 
progressive fracture has been added to the CEMCAN computer code. 
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The objectives of this paper are to describe the new features of the CEMCAN computer code and demonstrate its 
capabilities by simulating room-temperature tensile stress-strain curves up to failure of both unidirectional and angle-ply 
laminates of SiC/RBSN composites at room temperature. 


CEMCAN COMPUTER CODE 

In-house research at NASA Lewis Research Center in Cleveland, Ohio over the past two decades has resulted in 
several computer codes for composite micromechanics, mesomechanics and macromechanics. The primary intention of this 
research is to develop composite mechanics theories and analysis methods that range in scale from micromechanics to 
global structural analysis in one integrated code. The micromechanics theories are represented by simplified equations. 
Among the codes that have resulted from this effort are ICAN (Integrated Composite Analyzer) (ref. 1) and METCAN 
(Metal Matrix Composite Analyzer) (ref. 2). The ICAN computer code which is developed for the analysis of polymer 
matrix composites uses a representative volume element (RVE) that consists of fiber and matrix. A perfect bonding 
condition is assumed between the two constituents. The unit cell or RVE for METCAN, which is developed for the 
analysis of metal matrix composites, consists of fiber, matrix and an optional third phase called interphase. The interphase 
is treated as a distinct region with distinct mechanical and thermal properties which can be different than those of fiber or 
matrix. Thus, the interphase can represent a reaction zone formed due to the chemical reaction between fiber and matrix 
or a fiber coating or a layer provided intentionally to prevent such a reaction. This approach has been taken one step 
further for the analysis of ceramic matrix composites. Although, the same square fiber array pattern is used as in the 
previous codes, the fiber in the unit cell or RVE is divided into several slices. The composite micromechanics equations 
are derived for the slices rather the RVE. The equivalent properties of the RVE are derived by the application of classical 
iflminfltp. theory. The advantages of fiber substructuring technique and related micromechanics equations were described in 
reference 3. These equations have been programmed into the computer code CEMCAN. The flow-chart of the CEMCAN 
computer code is shown in figure 1. The left side of the chart shows the various steps involved in the synthesis of the 
properties from slice to laminate level, while the right side of the chart shows the various steps involved in the decomposi- 
tion of the response from laminate to ply, to slice and then to microstresses. Composite micromechanics equations are used 
to determine the equivalent elastic properties of a composite material in terms of the elastic properties of its constituent 
material g The properties of interest are composite moduli, Poisson’s ratios, thermal expansion coefficients, thermal 
conductivities, heat capacities etc. The ply/fiber substructuring is shown schematically in figure 2. This technique allows a 
more accurate micromechanical representation of the composite material and at the same time, provides a greater detail in 

local stress distribution within a ply. ... . , 

Recently, a nonlinear analysis capability that incorporates local stress redistribution due to progressive fracture has 
been added to the CEMCAN computer code. The material nonlinearity is accounted for at the constituent (fiber/matrix/ 
interphase) level. For example, the constituent material properties as a function of temperature during the fabrication/ 
cool-down process can be easily accounted for while computing the residual stresses. This capability allows for the 
inclusion of residual stresses that arise due to the temperature difference between the processing and the use temperatures 
and the migmairh in the coefficients of thermal expansion of the fiber and matrix materials prior to the application of any 
loading and simulates the stress/strain state in the as-fabricated material. The other capability that has been added is 
microstress redistribution as the fracture initiates and propagates in a laminate. If a constituent fails in a slice, it unloads 
and its microstresses are redistributed to the unfailed regions. This capability is described in more detail in a later section. 


MULTIFACTOR INTERACTION RELATIONSHIP 

Research over the last two decades at NASA Lewis Research Center to simulate the behavior of composite 
ma terial g has culminated in a unified law to describe the material degradation behavior. This unified law, referred to as the 
multifactor interaction relationship (MFIR), takes into account the degradation effect on constituent material behavior due to 
environmental, fabrication and load effects. A generic form of the above mentioned MFIR is given by. 
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where 

P current property 

P 0 reference property 

V current value of the primitive variable representing load/environmental effect 
superscript 

a exponent for a given effect 

subscript 

f final value of the variable 

0 reference value of the variable (e.g. for temperature effect, V 0 can be room-temperature) 

Any number of effects and their interaction on a given material property can be described. The exponents are 
determined from available experimental data or sometimes estimated from anticipated material behavior based on an 
educated guess. The nonlinearity in the overall composite behavior also arises due to micro fracture and subsequent 
microstress redistribution as the fracture initiates and propagates. In the present work, the constituent properties namely the 
modulus and coefficient of thermal expansion are considered a function of temperature only. Hence, for the modulus for 
example, the above equation (1) reduces to: 



( 2 ) 


where 

E is the modulus at temperature T 

Eo is the reference modulus at temp T 0 (usually modulus at room-temperature) 

T f is the final temperature where the modulus is zero 

MICROSTRESS REDISTRIBUTION DUE TO PROGRESSIVE FRACTURE 

A unique capability to carry-out a progressive fracture analysis with microstress redistribution resulting from 
progressive fracture has been added to the CEMCAN computer code. Once the microstresses are computed in all the 
constituents of a slice, they are compared with the current corresponding strengths. The maximum stress criterion for 
failure is used in the current analysis. The constituents are assumed to be void or flaw free. Although, for brittle 
materials, the strength values may depend upon the size and distribution of flaws, this failure criterion does not account for 
such effects. A statistical/probabilistic analysis should be carried out to accurately account for the effect of flaw size and 
distribution on the strength of a constituent. If a constituent stress has exceeded the corresponding strength value, then the 
constituent is assumed to have failed and the corresponding modulus of that constituent is reduced to an almost negligible 
value. A slice is assumed to have faded in the longitudinal direction if there is a fiber failure in that slice in the 
longitudinal direction. A matrix failure in the transverse or shear direction means the slice has failed in the transverse or 
shear direction respectively. Once a slice fails, all the microstresses that are carried by that slice up to that load step are 
redistributed in the slices that have not failed. The load redistribution due to progressive fracture is shown schematically in 
figure 3. The load that was carried by the failed slice is then added appropriately to the laminate load (as shown in fig. 3), 
and a laminate analysis is carried out again. This process is continued for a given load step until an equilibrium 
(convergence) state is reached between the applied load and a damage state. Presently, the convergence is checked for 
laminate mid-plane strains, ply and slice strains. Once all the strains are within 5 percent of the previous iteration, 
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convergence is assumed to have been reached. This mechanism is mainly responsible for the overall nonlinear behavior in 
the laminate tensile stress-strain curves at room-temperature. 


RESULTS/DISCUSSION 

The composite system used here is SiC/RBSN (silicon carbide SCS-6 fibers in a reaction bonded silicon nitride 
Si 3 N 4 matrix) as mentioned before. The fiber and matrix mechanical and thermal properties, shown here in table I, have * 

been provided in reference 4. This composite system is known to have a very compliant and weak interphase. The 
thickness of the interphase, a combination of outer carbon layer of the fiber and any reaction zone, is estimated to be 
approximately 3 percent of the fiber diameter (0.03 x 142 = 4.3 pm). However, the thermal and mechanical properties of 
the interphase are not known. The nominal fiber volume ratio for this composite was estimated in the range of 0.3 to 0.34 
in reference 5. Since it is well known that the interphase properties have very little effect on the longitudinal modulus of a 
unidirectional composite (ref. 6), the nominal fiber volume ratio (fvr) is calibrated to be 0.36 by comparing the measured 
and predicted values of longitudinal modulus for a unidirectional composite (fiber and matrix modulus are known). It 
should be mentioned here that in the simulation process, the formation of the interphase degrades the fiber. Hence, the 
effective fiber diameter is reduced and so is the fiber volume ratio but the matrix volume ratio remains unchanged. For 
example, in the present analysis, the gross fiber volume ratio is 0.36, fiber diameter is 142 pm, interphase thickness is 
3 percent of the fiber diameter, then effectively fiber volume ratio reduces to 0.318, interphase volume ratio is 0.042 and 
matrix volume ratio remains unchanged at 0.64 (1 to 0.36). Once the nominal fvr is known, the interphase modulus is 
calibrated by comparing the measured and predicted values of transverse modulus for a unidirectional composite since the 
transverse composite modulus is quite sensitive to interfacial conditions. Alternatively, the measured value of in-plane 
shear modulus for a unidire ctional composite can also be used to calibrate the interphase shear modulus. Once calibrated, 
these values were used to predict some mechanical properties for two laminates [0 2 /90 2 ] s and [±45 2 ] s and these predictions 
were compared with experimentally measured values. Results are shown in table n. The comparison is excellent, except 
for V] 2 for [±45 2 ] s laminate and all the predictions are within the scatter of the measured values. This calibration process 
is explained in little more detail in reference 4. A linear laminate analysis was also performed on a unidirectional 
composite and microstresses were computed for different types of applied loading. Although the results are not shown 
here, those resulting microstresses were compared with the predictions of a detailed three-dimensional finite element 
analyses and the results were in very good agreement. 

In the present work, the tensile stress-strain behavior up to failure for unidirectional on- and off-axis specimens 
and two angle-ply laminates is predicted and compared with experimental stress-strain curves at room temperature. The 
fabrication of these composites is explained in Ref. 5. These composites are processed (i.e., nitridation is performed) 
between 1200 to 1400 °C. However, after performing several analyses, the stress free temperature in the present 
simulations is assumed to be 900 °C in order to obtain the observed failure modes in these composites. This suggests that 
at higher temperatures (900 to 1200 °C) some metallurgical mechanism is operative that prevents stress build-up above the 
assumed stress-free temperature (900 °C). Such observations have been made before with regards to other ceramic matrix 
composites (ref. 7). In the present simulations, the composites are cooled down to room-temperature (20 from 900 °C) 
prior to application of any loading. It is also mentioned that the maximum use temperature of these composites will be 
around 1200 °C and there is very little degradation (maximum 10 percent, ref. 8) in constituent modulus etc. even at the 
maximum use temperature compared to the reference property at room temperature. To take that into account, the 
exponent for the temperature dependence of the properties (modulus and coefficient of thermal expansion) is taken as 0.25 
for the fiber and 0.1 for the matrix, final temperature (T f , equation (2)) is taken as 2500 °C for the fiber and 2200 °C for 
the matrix material. These assumed values would correspond to the observed degradation in the constituent material 
properties, namely modulus and thermal expansion coefficient for the range of the use temperature for these material. 

Room-temperature stress-strain curves for both [0] 8 and [90] 8 composites are shown in figure 4. Based on the 
observed tensile strengths of [0] 8 and [90] 8 composites, matrix in-situ tensile strength was calibrated as 100 MPa (14.5 ksi) 
and average fiber bundle strength as 2 GPa (285 ksi) by comparing the knee in the stress-strain curve of [0] 8 composite 
i.e„ the deviation in the linear behavior and the fiber strength is calibrated by comparing the final failure stress of the same 
composite. These values agree well with the fiber and matrix tensile strengths reported in reference 5. Once calibrated, 
the same strengths values were used in subsequent simulations. For the [0] 8 composite, stress-strain curve is linear initially 
as both fiber and matrix are intact Upon further loading, the matrix starts cracking until saturation is reached. Conse- 
quently, all the load is carried by the fibers only, which is represented by the second linear region of the stress-strain curve. 

Final fracture occurs when the fibers fail in all the plies. Such an observation agrees very well with the optical 
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observations marie of specimens loaded to different stress levels (ref. 5). The [90]g composite shows linear behavior until 
failure. It failed with the failure of matrix in the transverse tension, at a relatively low stress level of 45 MPa. 

In contrast, the off-axis laminates [10] 8 and [45] 8 stress-strain behavior is largely linear up to failure as shown in 
figure 5. Such a behavior is characteristic of a brittle failure. The final fracture in these laminates is controlled by a 
combination of matrix/interphase normal and shear failure modes. Optical micrographs shown in reference 5 are in 
agreement with such failure modes. Room-temperature stress-strain curve to failure for a [±45] s laminate is shown in 
figure 6. The shear strength of the matrix is estimated to be half of the matrix tensile strength to match the ultimate stress 
of [±45 2 ] s laminate. This laminate also fails in a combination of shear and normal failure modes and shows the charac- 
teristics of a graceful failure. Once again, the predicted strength of this laminate matches well with the experimentally 
measured value, and the optical micrographs in reference 5 confirm the above mentioned failure modes. The unloading 
part of the stress-strain curve for a [±45 2 ] s laminate is not predicted. 

The room-temperature stress-strain curve to failure for a cross-ply laminate [0 2 /90 2 l s is shown in figure 7. Up to a 
stress level of 150 to 200 MPa, the curve is linear and shows no damage. The initial modulus agrees well with the experi- 
mental measurement. Damage initiates, first in the form of fiber matrix debonding in the 90 plies and subsequently as 
matrix cracking in the 0° plies. In the second linear region of the stress-strain curve, the 90° plies are fully cracked and do 
not contribute to the stiffness in the loading direction and the matrix is fully saturated with cracks in the 0° plies and does 
not contribute to the stiffness in the loading direction either. All the load at this point is earned by fibers only in the 0 
plies. However, the slope of the second linear portion of the stress-strain curve is higher than that of the observed, stress- 
strain behavior. The experimental curve suggests that all the fibers in the 0° plies are not contributing to the stiffness. One 
possible explanation is that there is probable damage due to processing in the fibers of 0° plies. All the 0° fibers are not 
contributing to the stiffness in the second linear portion of the stress-strain curve. Comparing the slope of the second linear 
portion of the stress-strain curve, it would appear that the fiber volume ratio of 0° plies is approximately half, i.e., experi- 
mental stress-strain curve slope of the second linear portion is approximately half of the slope of the predicted curve. The 
fiber lay-out for this laminate is shown in figure 8. It is possible that during the processing/composite fabrication, the 
application of high pressure at relatively high temperature, will result in impingement of 90° fibers on to the inner 0° 
fibers, resulting in some damage to the inner 0° fibers. If both the inner fibers were damaged, that will be two fibers out 
of four in the loading direction, and would corroborate with the earlier observation that only half of the fibers in loading 
direction seem to be effective. Damage to the fibers due to processing could result in a shift in the mean strength of fibers 
(fig. 8) causing premature fiber failure. Although the predicted stress-strain curve accounts for the residual stresses 
resulting from the processing of the composite, it assumes no damage in the 0° fibers and thus results in a much stiffer 
behavior and a higher strength value. In fact, the experimentally measured strength value for this laminate (290 MPa) is 
bound by the predicted strength with no damage in the fibers, 390 MPa (upper bound) and with both inner 0° fibers 
damaged, 195 MPa (lower bound). 


SUMMARY/CONCLUSIONS 

Recent upgrades to the CEMCAN computer code including nonlinear analysis capability with microstress 
redistribution due to progressive fracture are described. Collectively, the results presented show that CEMCAN provides 
sufficient flexibility to the user to study various aspects in the deformation/fracture behavior of ceramic matrix composites. 
Based on the results presented and discussed, it can be concluded that: 

1 . Unidirectional [0] 8 composite loaded along longitudinal direction displays linear stress- strain behavior up to first 
matrix cracking stress. Beyond this point, matrix is ineffective and all the load is carried by the fibers until fracture. 

2. Off-axis and angle-ply laminates tensile strength is governed by a combination of matrix/shear failure. 

3. There is probable fiber damage/breakage due to processing in the cross-ply [O^C^ls laminate. 

4. CEMCAN computer code provides the user with a “virtual test laboratory,” in that the new candidate materials, 
different laminate configurations or service environments can be evaluated in a fraction of time and cost, with only a 
minimum number of tests needed for the calibration of certain material properties for a given material system. 
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TABLE I.— CONSTITUENT PROPERTIES OF SIC/RBSN AND PROPERTIES OF [0] g 
COMPOSITE CONSTITUENT PROPERTIES 


Property 

SiC (SCS-6) fiber 

RBSN matrix 

Interphase 

Modulus, GPa (Msi) 

390 

(56.6) 

no 

(15.95) 

3.5 (calculated) (0.5) 

Poisson’s ratio 

0.17 


0.22 


0.22 


Shear modulus, GPa (Msi) 

117 

(17) 

45 

(6.5) 

1.4 

(0.2) 

Coefficient of thermal 
expansion. l(rVc (lO^F) 

4.1 

(2.3) 

2.2 

0-2) 

Z0 

(i.i) 

Thermal conductivity, 
W/m-K (Btu/ft-h °F) 

22 

(12.7) 

5 

(2.9) 

ZO 

(1*2) 


UNIDIRECTIONAL SIC/RBSN [0] g COMPOSITE FIBER 
VOLUME RATIO 0.36* 


Property 

Predicted 

Measured^ 

Bn, GPa (Msi) 

189 

(27.4) 

193±7 

(28±1) 

E^, GPa (Msi) 

68.3 

(9.9) 

69±3 

(10±0.4) 

Gj 2 » GPa (Msi) 

27 

(3.9) 

31±3 

(4.5±0.4) 

V 12 

0,21 


0.2 



a Note: The measured value of longitudinal modulus for a 
unidirectional composite, Ej j, was used to compute fiber 
volume ratio (fvr), while the transverse modulus, E^, was 
used to compute the interphase modulus. 

Reference 3. 
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TABLE n — COMPARISON OF PREDICTED (CEMCAN) AND MEASURED 
ELASTIC PROPERTIES OF SIC/RBSN LAMINATES 


Property 

[0 2 A*0 2 ] S lami 11 ^ 

[±45 2 ] $ laminate 

Predicted 

Measured 3 

Predicted 

Measured 1 

E n , GPa (Msi) 

129 (18.7) 

124±6 (18±0.9) 

78.6 (11.4) 

78±3 (11.3±0.4) 

& 22 , (Msi) 

129 (18.7) 

124±6 (18±0.9) 

78.6 (11.4) 

78±3 (11.3±0.4) 

G| 2 » GPa (Msi) 

27 (3.9) 

31 ±2 (4.5±0.3) 

58 (8.4) 


V I2 

0.112 

0.12 

0.46 

0.36 


■=> as ■=> 



Slice 


Figure 1 . — Integrated analysis approach embedded in CEMCAN 
computer code. 
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Horizontal slicing Vertical slicing 

(along the 2-2 direction) (along the 3-3 direction) 

(Note: 1-1 is along the fiber direction) 



fiber, matrix, and interphase 

Figure 2. — Ply/fiber substructuring concepts for ceramic matrix 
composites micromechanics. 



• Failure based on maximum 
stress criterion 

• Fiber failure governs the 
failure of slice in 

the longitudinal direction 

• N x * * N x + s$iice 11 x 
(slice thickness) 



• Matrix/interphase failure 
governs the slice failure in 
transverse (shear) loading 

• Ny' = Ny + Ssijce 22 x 

slice thickness) 


Figure 3. — Local stress redistribution due to progressive fracture. 
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Figure 4. — Room -temperature stress-strain curves 
to failure of [0] 8 and [90] 8 SiC/RBSN composite. 



Figure 5. — Room -temperature stress-strain curves to 
failure of [1 0] 8 and [45] 8 laminates of SiC/RBSN, 



Figure 7. — Room-temperature tensile stress-strain curve to 
failure of [Oa/SO^s cross-ply laminate of SiC/RBSN. 




Loading direction 



Figure 6. — Room-temperature stress-strain curves to failure of 
[+45/~45] s SiC/RBSN laminate. 
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► Processing can induce 
damage in the fibers 

and shift the mean strength 

► Experimentally measured 
value of [Og/SO^s laminate 
strength is bound by no 
fiber damage and fiber 
damage in two [0] plies 


Fiber strength 


Figure 8. — Fiber layout for [(>2/90 Js laminate. 
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